Theory of superconducting pairing near the mobility edge 
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We develop a theory of a pseudogap state appearing near the superconductor-insulator transition 
in strongly disordered metals with attractive interaction. We show that such an interaction combined 
with the fractal nature of the single particle wave functions near the mobility edge leads to an 
anomalously large single particle gap in the superconducting state near SI transition that persists 
and even increases in the insulating state long after the superconductivity is destroyed. We give 
analytic expressions for the value of the pseudogap in terms of the inverse participation ratio of the 
corresponding localization problem. 



A rapidly growing number of experimentsi^^^^i^ 
on various disordered superconductors show that a novel 
phase appears in the vicinity of the superconductor- 
insulator transition as sketched in the low temperature 
phase diagram in Fig.fJ On the superconducting side of 
the transition a relatively small magnetic field suppresses 
the superconductivity leading to an insulator character- 
ized by a large thermally assisted resistance with a hard 
gap. Upon a further increase of the magnetic field the 
resistance and the gap drop 4 i 5 i 6 i 7 i 8 . It is tempting to 
explain these data by a formation of localized Cooper 
pairsiS - in this picture the superconductivity is due to 
a fragile coherence between localized Cooper pairs, while 
the energy of single electron excitations is much larger 
and remains finite even when the coherence (and thus 
the superconductivity) is destroyed. Experimentally, this 
behavior is observed only for a limited range of disorder 
strengths, a± < a < era, characterized by the conductiv- 
ity a at room temperature. At weaker disorder, a > 
the destruction of superconductivity by a magnetic field 
leads, as usual, to a formation of a normal metallic state 
without any noticeable localization effects. At larger dis- 
order a < a 1 superconductivity is replaced by a gapped 
insulator even at B — 0. Usual insulator with a variable- 
range hopping conductance is recovered at larger disor- 
der, a < (To. 

The hypothesis of preformed Cooper pairs is further 
confirmed by the behavior of disordered superconductors 
at higher temperatures. On the insulating side of the 
transition in thick (effectively 3D) films one observes^ 
Arrenius behavior, R(T) oc exp(T//T), at low tempera- 
tures. The experimental value of the activation energy, 
Tj, is somewhat larger than the superconducting gap in 
less disordered samples and grows with the disorder. For 
instance in InO^ films Tj = 2 — 15 K, while the criti- 
cal temperature T c — 3.3 A' in best samples 1 ^. How- 
ever, at higher temperatures (T > 10K for InOx films 
studied in Ref. |3J) this behavior is replaced by Mott's 
variable range hopping R(T) ~ exp(T M /T) 1 / 4 . Thus, 
qualitatively the temperature plays the same role as the 
magnetic field, suppressing an unconventional insulator 
formed near superconductor-insulator transition. This 



can be easily understood if the insulating pseudogap is 
due to preformed Cooper pairs with a pairing energy 
somewhat larger than the superconductive gap in a less 
disordered samples. 
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FIG. 1: Schematics of the phase diagram of disordered su- 
perconductors in the vicinity of the superconductor-insulator 
transition as a function of the magnetic field, B, and disorder 
strength that is characterized by room temperature conduc- 
tivity, a ~ UfI- In the shaded region the resistiviy exhibits 
activated behavior with a well defined gap Tj. 

The purpose of this Letter is to provide a justifica- 
tion for this picture and to construct the theory of local- 
ized Cooper pairs starting from a semi-microscopic model 
that contains only low energy electrons with BCS-type 
attraction and a strong random potential which leads to 
Anderson localization of single-particle states. An im- 
portant ingredient of the problem is the fractal structure 
of single-particle wave functions, which sets in on a mi- 
croscopic scale Lo determined by the disorder potential 
and extends to the localization length L\ oc . 

It is well established that disorder has little effect 
on conventional superconductivityiiiiS in the absence of 
Coulomb repulsion. On the mean-field level (neglect- 
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ing thermal fluctuations), T c and the self-consistent BCS 
gap A are affected by the localization and start to de- 
crease only when the effective level spacing in the lo- 
calization volume, Sl = l/(z/o£? oc ), becomes compara- 
blo 11 T 13 i 14 i 15 i 16 to A (here v is the density of orbital 
states at the Fermi level). Intuitively, it is clear that the 
thermodynamics of a strongly disordered superconductor 
built on localized electron wave functions is practically 
the same as that of a system separated into compart- 
ments (grains) of size L\ oc . In each grain Sl serves as 
a low energy cut-off; when it exceeds A the supercon- 
ductivity is destroyed. However, the attraction between 
electrons persists until Sl exceeds a high energy cut-off 
which is typically the Debye frequency, cod- The attrac- 
tion results in pairing of electrons with opposite spins 
occupying the same localized state. Assuming that the 
pairing interaction is local and is characterized by a di- 
mensionless coupling A, the pairing energy, 2Ai, is pro- 
portional to the inverse volume occupied by the state, 
Ai ~ A/(^oKff). It is crucial for the following that this 
volume is paramctrically smaller than the 'localization 
volume' if oc due to the fractal nature of the wave func- 
tions near t he mo bility edge V e g ~ Lq(L\ oc / L ) D2 , where 
D2 < 3, cf. Il7ll8( Then the energy scale Ai is paramet- 
rically larger than Sl- Thus, the strength of the disorder 
is characterized by two relevant energy scales, Sl and Ai. 
Here we focus on two limiting cases 

A < S L <C oJd (I) 
5 L < A«Ai (2) 

As shown below in the first regime the single-particle 
gap remains large leading to a gapped insulator. In the 
second regime the superconductivity persists, but the 
single-particle gap Ai is much larger than in conven- 
tional superconductors. In particular, the ratio of this 
gap to T c can be anomalously large. Moreover, in the 
latter regime the superconducting phase transition is pre- 
empted (upon T decrease) by a formation of a "pseudo- 
gaped" insulating state similar to the one observed in 
underdoped cupratesi^. 

We assume that a BCS-type attraction between elec- 
trons at high energy scales ~ ujl> is only weakly affected 
by localization of electron wave functions provided that 
Sl -C wi). The theoretical reason for this is the scal- 
ing argument that says that physics at short scales and 
large energies is not affected by localization as long as 
kFLi oc 3> 1. This is further collaborated by an experi- 
mental observation^ that T c changes little until one gets 
into a deeply localized regime. In a fermion system with 
weak attraction one can leave only the pair interaction 
terms in the Hamiltonian leading to the usual BCS model 
in the basis of localized electron states^: 

H = € 3 c ]a c 3<y ~ — Yl X 3kM jk c^c f n c k1 c kl (3) 

where €j is the single-particle energy of the state j and the 
matrix elements Mjk = J drip^(r)?pl.(r). In Hamiltonian 



© we will also distinguish the off-diagonal dimensionless 
coupling constant Xjk = A for j =/= k from the diagonal 
one Xjj = A for the reasons detailed below. 

To study the model © we need the statistics of single 
electron states near the 3D Anderson mobility edge. We 
expect that on average off-diagonal matrix elements, Mjk 
with j k have the usual scaling l/£f oc for the states 
localized at distances smaller than L\ oc from each other. 
In other words we assume that correlations Mjk be- 
tween wave functions of different states are much smaller 
than individual "inverse participation ratios" (IPRs), 
Mjj = J drtpj (r) , which have the meaning of the inverse 
volumes of localized states, V e s- In fact, extensive nu- 
merical studies^ii^ show that a typical IPR, M, scales 
as M w L 3 (L loc /L )- D2 , with D 2 = 1.30 ±0.05. 

As discussed above, the thermodynamics of this model 
should be similar to a collection of grains of size L\ oc with 
Hamiltonians 

H = J2 e i c ]<y c i° ~ 9]kc)^c] iCk] Cki (4) 

where the average level spacing is (e/+i — tj) = Sl, the 
off-diagonal part of the interaction has the usual form 
for a grain, gjk = g = XSl for j ^ k, while the diag- 
onal part has an extra term G so that gjj = g + G = 
X(L\ oc /Lq)~ D2 /(vqLq) due to the fractal nature of the 
wavefunctions. We have checked that the results de- 
scribed in this Letter are the same (up to numerical fac- 
tors of order one) for models and J3J. 

We begin with the insulating region Q where Cooper 
interaction can be treated perturbatively^S. In the first 
order perturbation theory the minimum energy (counted 
from the chemical potential), Ai, required to add a single 
electron to the system is 

A / \ 

1 = 2v Ll(L loc /L ) D2 (5) 

The coupling constant A appearing in this equation is 
not renormalized by Cooper loops: indeed, the Hamilto- 
nians © and Q do not mix unoccupied and doubly oc- 
cupied (unblocked) states with singly occupied (blocked) 
ones. Blocked states do not participate in pair scatter- 
ing and contribute only their single-particle energy to the 
Hamiltonian. The additional contribution coming from 
the diagonal term in @ is GN pa i rs , where N pa i rs is the 
total number of pairs. Using the fact that the total num- 
ber of electrons is 2N pairs + Nb , where Nb is the num- 
ber of blocked states and shifting single-particle levels 
ej — » Cj + G/2, we can rewrite Hamiltonian (0} as 

ja j.k B ^ ' 

(6) 

Here the first two terms act only on unblocked states, 
while the summation in the last term is over decoupled 
blocked states. The renormalization affects only the in- 
teracting part which is of a standard BCS form. This 
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statement is valid for the reduced Hamiltonians (USUI) 
which contain pair-wise eigenstates only. In general, one 
would say that our phenomenological Hamiltonians (|3l4jl 
contain two different coupling constants originating from 
the competition between Coulomb and electron-phonon 
interactions - standard BCS coupling A and the addi- 
tional "diagonal" one, A, which characterizes interaction 
of two electrons in the same localized state. 

Applying the same arguments to the Hamiltonian J3jl, 
we get a single-particle gap that varies from one level 
to another: Aj = \Mjj / (2vq). The average density of 
single electron states in a large sample is controlled by 
the IPR distribution V{M) 



2ev /X 



V(M)dM 



(J) 



Scaling theory of localization predicts that near the mo- 
bility edge V(M) acquires a scale- invariant form and 
this is indeed what was observed in numerical studies 1 ^. 
Moreover, these data indicate that the distribution V(M) 
decreases fast at M/M — > so that extended states occur 
too rare to smear up a gap-like behavior in the density 
of states (J7J with the gap value given by Eq. JSJ). 

We emphasize that the DoS Q does not contain any 
"coherence peak" above the gap. The same qualita- 
tive behavior of the DoS in the insulating region was 
obtained 24 by solving numerically BCS mean field equa- 
tions for the disordered "negative-U" Hubbard model. 
Evidently, the spectral gap Ai should be associated with 
the experimentally determined 1 ^ activation energy Tj. 
The gap dependence predicts a moderate increase 
(oc l/ii c) of Tj with disorder strength in a complete 
agreement with the experimental datap», see Fig. 2. Note 
that attempts to fit the data with a gap that scales as 
8l k 1 /Lf oc have failed spectacularly 1 . 




(J~(k F l)(au) 

FIG. 2: Experimental values of the gap, Tj (boxes) and a fit 
to the equation © with L; oc ~ (a — <t c ) _1 that involved only 
the overall scale adjustment^. 

We now turn to the parameter region (J2J where one ex- 
pects a global superconductive coherence. Namely, the 
condition Sl -C A assures that a large number of eigen- 
states t/jj with energies down to ej ~ A overlap signif- 
icantly in the real space, so one can use the standard 
BCS-type self-consistent gap equation as first argued in 
Ref . M 



However, due to an extra cost of creating an unpaired 
electron, the energy required to add an electron in this 
regime is no longer equal to the self-consistent gap A as 
in the usual BCS theory. Indeed, because the addition of 
an electron increases the number of blocked states in (JBJ 
by one, we have 



A 1 = A + - 



A 



2v L 3 (L loc /L ) D * 



(8) 



In order to evaluate the energy of pair excitations, i.e. 
the excitations of the condensate with no unpaired elec- 
trons, it is convenient to rewrite Hamiltonian (0 in terms 
of Anderson's pseudospin-1/2 operators 21 - defined on the 



subspace of unblocked states as 
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1. We have 



ff = 5> 
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(9) 



In the BCS mean-field approximation eigenstates of © 
correspond to all spins being parallel or antiparallel to 
the effective magnetic field Bj = — (2A,0, 2ej), where 
A = gJ2j( s j): leading to a self-consistency condition 



9 , 



A 2 



(10) 



j — 1 if the spin is against the magnetic field and 
1 otherwise. 



where e 

In the ground state all spins are along their mag- 
netic fields, ej — 1 for all j. The lowest pair excitation 
corresponds 21 * 2 ^ to a single spin flip at the Fermi level 
eo = and has an energy A2 = 2A. It also modifies 
the gap equation since now one of ej becomes negative, 
eo = — 1. Because the number of states in the localiza- 
tion volume Lf oc is finite, this modification results22iS in 
a negative correction to the gap: i5A = —5^/2. Thus 



2A 



2A 



(11) 



loc 



here Ao is the T = gap at large L/ oc , numerical coeffi- 
cient c = 1 for model and c ~ 1 for model @. Note 
that pair excitations with gap A2 are spinless, whereas 
single-particle excitations (gap Ai) carry spin 1/2. Thus 
Ai is the spin gap, while A2 is the energy gap. 

In the regime the typical single-particle gap Ai 
remains much larger than the gap for pair excitations A 2 . 
This leads to the exclusion of single-particle excitations 
from the thermodynamics at T Ai, which results in a 
modified mean-field result T C /A(T = 0) = 2[T C /A(T = 
0)]bcs- The actual dependence of T c on the localization 
length is not easy to find due to enhanced role of thermal 
fluctuations near the mobility edgei^iii^. 

We emphasize that A(T = 0) does not characterize 
the spectral gap in the present situation. The most nat- 
urally measured (e.g. by tunneling conductance or opti- 
cal conductivity experiments) spectral parameter is the 
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single-particle gap Ai, which according to Eqs. i|8l2[l is 
much larger than T c . An anomalously large ratio Ai/T c 
leads to the insulating trend of the resistivity versus tem- 
perature behavior in the intermediate temperature range 
T c < T < Ai. This was observed in many strongly dis- 
ordered superconductors and is especially well-known for 
underdoped cuprates^. 

The quantitative similarity between R(T, B) behav- 
ior in InOa; films and underdoped cuprates noticed in 
Ref [2(1 (see also Ref. makes us believe the pseudogap 
in underdoped cuprates may have a similar origin: pair- 
ing of electrons on localized states. The important differ- 
ence of the cuprates is the d-wave symmetry of the pair- 
ing. We plan to address this issue in future work. Sup- 
pression of single-particle density of states in the pseudo- 
gap regime can be observed by measuring the differential 
tunneling conductance which we expect to follow the uni- 
versal IPR distribution V(M) studied in Ref. 0: 

^^Lv. ^' 2 *"^ (i2) 

As emphasized above, the generic features indicat- 
ing a formation of the novel insulating state near 
superconductor-insulator transition (such as the resistiv- 
ity maximum as a function of magnetic field or disorder) 
were observed in many different materials. However, the 
details may vary: in InO^ thick (essentially 3D) films the 
resistivity above the transition has an activated behav- 
ior while in very thin (2D) films2*2& it follows an Efros- 
Shkovskii law. The theory developed in this Letter is 
applicable only to 3D systems but it seems likely that a 
similar physics leads to a weaker but observable effects 
in 2D films. 



In conclusion, weak Anderson insulators with Cooper 
attraction are shown to possess hard insulating gap 
whose magnitude is determined by the IPR statistics near 
the mobility edge. Although this gap is of a supercon- 
ducting origin, it does not lead to a coherence peak. In 
the ground-state of this insulator all electrons are paired 
on individual localized eigenfunctions. When the Fermi- 
level approaches the mobility edge, superconductive cor- 
relations develop between localized pairs. Key features 
of the predicted superconductive ground-state are an un- 
usually large (compared to T c ) single-particle excitation 
gap (spin gap) and a pseudogaped regime at tempera- 
tures about T c . These unusual features stem from the 
fractal nature of localized eigenstates near the mobility 
edge. We assumed that V^( r ) of the relevant nearly- 
critical eigenfunctions have weak mutual correlations. If 
the correlations are in fact strong^, the description of the 
insulating phase is not affected but the ratio Ai/A 3> 1 
will be significantly reduced. Finally, our theory is based 
on a phenomenological assumption of an attractive sign 
of the diagonal coupling constant A > 0. Another inter- 
esting situation is realized when A > 0, but A < 0. In this 
case we expect a formation of "effective magnetic impu- 
rities" and strong suppression of T c . We plan to consider 
this situation separately. 
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